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We study the heat transfer between two finite quantum systems initially at different 
temperatures. We find that a recently proposed fluctuation theorem for heat exchange, 
namely the exchange fluctuation theorem [C. Jarzynski and D. K. Wojcik, Phys. Rev. Lett. 
92, 230602 (2004)], does not generally hold in the presence of a finite heat transfer as in the 
original form proved for weak coupling. As the coupling is weakened, the deviation from the 
theorem and the heat transfer vanish in the same order of the coupling. We then discover 
a condition for the exchange fluctuation theorem to hold in the presence of a finite heat 
transfer, namely the commutable-coupling condition. 

We explicitly calculate the deviation from the exchange fluctuation theorem as well as 
the heat transfer for simple models. We confirm for the models that the deviation indeed 
has a finite value as far as the coupling between the two systems is finite except for the 
special point of the commutable-coupling condition. We also confirm analytically that the 
commutable-coupling condition indeed lets the exchange fluctuation theorem hold exactly 
under a finite heat transfer. 



The development of the modern techniques of microscopic manipulation enables 
us to treat small systems, for example, nano-devices and molecular motors. In such 
small systems, classical thermodynamics is not well applicable to quantification of 
the heat flow or the work. At the nano- and micro-scales, the available thermal 
energy per degree of freedom is comparable to the energy of the small systems. This 
thermal energy enhances fluctuations, whose effects measurably appear in such small 
systems. Thus, we cannot apply classical thermodynamics to these small systems. 
We need a substitute framework if we try to design or control nano-devices and 
molecular motors as macroscopic heat engines. 

Since the pioneering work by Evans et al., 1 ^ researchers have developed various 
results collectively known as the fluctuation theorem. The fluctuation theorem is 
roughly summarized as an equation which relates the probability p(f2) of observing 
the entropy increase i?, to the probability p(—f2) of observing the entropy decrease 
of the same magnitude: 



The definition of the entropy increase Q depends on the dynamics of the system 
under consideration. However, the fluctuation theorem has been established for 
several class of systems with reasonable definitions of Q. When the fluctuation 
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theorem is applied to the transient response of a system, the theorem is referred to 
as the transient fluctuation theorem. 2 ) 

For more than 15 years, many versions of the fluctuation theorem have been 
presented for a variety of classical-mechanical situations: thermostated systems, 1 ' 3 ) 
stochastic systems, 4 )' 5 ) and externally driven systems. 6 )~ 8 ) These classical- mechanical 
versions of the fluctuation theorem have been reviewed in Refs.9), 10). Some of these 
theorems were verified experimentally. n )~ 13 ) 

When the system becomes further smaller, quantum effects may become signif- 
icant. It is, however, not straightforward to extend the afore-mentioned classical- 
mechanical fluctuation theorem to quantum-mechanical systems. The crucial dif- 
ference of the quantum systems from the classical systems is an essential role of 
measurements. In order to generalize the fluctuation theorems to quantum systems, 
we need to identify the entropy, the work, and the heat that are measured in the 
quantum-mechanical context. 

There have been two attempts to do this: first, defining operators to represent 
the heat and the work; second, measuring the system and using the measurement 
outcomes to represent the heat and the work. In general, the former attempt has 
led to quantum corrections to the classical results. 14 )~ 17 ) In the latter attempt, on 
the other hand, several versions of the fluctuation theorem have been shown without 
quantum corrections. 18 )~ 22 ) Both the heat and the work are defined as the difference 
between the results of two measurements, a two-point quantity. We refer to this 
attempt as the two-point measurement. 

The exchange fluctuation theorem (XFT) 23 ) was proposed in the framework of 
the two-point measurement by C. Jarzynski and D. K. Wojcik. 23 ) The theorem 
concerns the statistics of heat exchange between two finite systems of Hamiltonian 
dynamics, initially prepared at different temperatures. Let (3a and (3b denote the 
inverse temperatures at which the systems A and B are prepared, respectively. The 
symmetry relation of the exchange fluctuation theorem is expressed with the prob- 
ability distribution p T (Q) of the net heat transfer Q as follows: 

Pr(Q)=e A ^ Pr (-Q), (1-2) 

where A(3 = (3b — (3a is the difference between the inverse temperatures and r is 
the time duration between the two measurements. They referred to this symmetry 
relation as the exchange fluctuation theorem because the exchanged energy Q be- 
tween the two systems is regarded as a heat transfer. Equation (|l-2p implies that the 
average of e~ A/3 ® over the ensemble of realizations for any time interval r is unity: 

(e~ Af3Q ) T = 1. (1-3) 

This integral form of equality is a direct consequence of the exchange fluctuation 
theorem, and thus we refer to this equality as the integral exchange fluctuation 
theorem (IXFT). 

Jarzynski and Wojcik derived the exchange fluctuation theorem in both classical 
and quantum systems. We will discuss the exchange fluctuation theorem only for 
quantum systems in the present paper. The situation in which the exchange fluc- 
tuation theorem was "proved" 23 - 1 is quite simple; the two finite quantum systems 
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are prepared in equilibrium at different temperatures, and then placed in thermal 
contact with one another. There is no work resource such as an external field nor an 
external force. In this situation, we simply identify the energy increase of each sys- 
tem as a heat flowing into the system. Equations (jl-2p and (jl-3p were then "proved" 
in the case where the coupling between the two quantum systems is weak. 
To summarize the present paper, we find the following: 

1. The exchange fluctuation theorem does not generally hold in the presence of a 
finite heat transfer as in the original form proposed for weak coupling. 

2. If the Hamiltonian that couples the two systems commutes with the total Hamil- 
tonian, the exchange fluctuation theorem becomes an exact relation for arbi- 
trary strength of coupling. We refer to this condition as the commutable- 
coupling condition. 

In short, the exchange fluctuation theorem holds not because of weak coupling but 
because of the commutable-coupling condition. 

The present paper is organized as follows. In £J2J we first introduce the joint 
probability in the two-time measurement and derive a symmetry relation about the 
joint probability. Then, we show deviation from the exchange fluctuation theorem. In 
£}3l we introduce a condition on the coupling Hamiltonian for the exchange fluctuation 
theorem to hold exactly, which we refer to as the commutable-coupling condition. In 
£01 we show an explicit form of the deviation for simple models and demonstrate that 
the deviation from the exchange fluctuation theorem vanishes under the commutable- 
coupling condition. Conclusions will be presented in £j5j 



2.1. Joint probability 

In this section, we define a procedure of the two-point measurement and the 
corresponding joint probability. Assuming the time-reversal invariance, we derive a 
symmetry relation of the joint probability. 

We consider two finite quantum systems given by the Hamiltonian 



where Ha and Hb are the Hamiltonians of the systems A and B, respectively. The 
second term of the right-hand side of Eq. (|2-ip . H c , is the coupling Hamiltonian 
which describes the connection between the two systems and 7 is a factor controlling 
the coupling strength between the two systems. In the present paper, we consider 
the case where H c does not have any extra degrees of freedom other than in the 
systems A and B. 

Let \m a ) and denote an eigenstate of H a and the corresponding eigenvalue 
(a = A,B), respectively. We refer to the product states \rriA) <8> |wib) as \m) for 
simplicity. The Hilbert space is spanned by the set of \m). 

Here, we introduce the measurement procedure to discuss the heat exchange 
between the two systems over a finite time duration < t < r. We consider a 



§2. Exchange fluctuation theorem in quantum mechanics 



H =H + 7 # c , 
H =H A ® 1 B + 1 A <g> H B 





4 



S. Akagawa and N. Hatano 



projection measurement (ideal measurement) only, and thus the state of the system 
is projected onto the corresponding eigenspace of Hq = Ha + Hb after the energy 
measurement. Using the eigenstate |m), we can write the projection operator onto 
the eigenspace of the Hamiltonian Hq with the eigenvalue E m = E A A + as 

n m = \m)(m\. (2-3) 

If the Hamiltonian Hq has degeneracy, we distinguish the degenerate eigenstates 
with a quantum number A m as \m, A m ). In this case, the projection operator onto 
the eigenspace of Hq with the eigenvalue E m is given by 

n m = ^2\m,\ m )(m,\ m \. (2-4) 

Am 

In the following discussion, we assume the eigenvalues of Hq to be nondegenerate for 
simplicity. Although Jarzynski and Wojcik did not mention the degeneracy of the 
Hamiltonians in their original paper, 23 ) the degeneracy indeed causes no corrections 
to their discussions. 

For time t < 0, each of the systems A and B is separately connected to a heat 
reservoir at the inverse temperatures (3a and (3b, respectively, for sufficiently long 
time to reach its equilibrium state. The reservoirs are removed just before t = 0, and 
hence the initial state of the total system is given by the following product state: 

Pmit=PA®pB, (2-5) 

Pa =—y~, (2-6) 

where Z a = Tre _/3a ^ Q is the partition function of the system a (a = A,B). At 
t = 0, we perform the first measurement of the energy of each system. Suppose that 
we obtained the outcome (E^ A , E^ B )- Then, the state of the systems is projected 
onto the corresponding eigenstate of Hq = Ha <8> 1_b + <8> Hb' \m) = \mA) <8> 
with the probability 

(m\ p in i t \m) = — — . (2-7) 

^A^B 

The coupling between the two systems is turned on just after the first mea- 
surement at t = 0. Then, the total system evolves according to the von Neumann 
equation from t = to t = r, where the time evolution is described by 

U(t) = e - i (^0+7^c)i/^. 

At t = t, we separate the two systems again and measure the energy of each system. 
Suppose that we obtained the outcome (E A A , E^ B ). Then the state of the systems 
is projected onto the corresponding eigenstate of Hq: \n) = \ua) <8> |ne). The joint 
probability with which we observe {E^ A , E^ B ) at time t = and {E A A , E^ B ) at 
time t = t in the above process is given by 

e -flAE* A ~P B EB B 

P(m, n|r) = — T„(r), (2-8) 

^A^B 
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where 



T„(r) = \(n\U(r)\m}\ 2 (2-9) 



is the transition probability from \m) to |n). Using the projection operator 77 m , we 
can write the joint probability (|2-8H as 



P(m,n\r) = Tr ^i7„C/(r)77 mAnit 77 m C/t(r) 

= TV (n n U(r)p init n m U^(T)) . (2-10) 



In the degenerate case, the joint probability is given by the same form as Eq. ([2-lOp . 

For derivation of the exchange fluctuation theorem, we hereafter assume the 
time-reversal invariance of the system. The time-reversal operation in quantum 
mechanics is described by an antilinear operator <9. 24 ) The time-reversal invariance 
of the system is then expressed as 

[O,H A ] = 0, [O,H B ]=0, (2-11) 
[e,H c } = 0. (2-12) 

The commutation relations in Eq. (|2-llj) are equivalent to the commutation relation 
between the projection operator 7T m and the time-reversal operator 0: 

[O,n m ]=0. (2-13) 

By inserting the identity 1 = O~ l between U(t) and /9; n it in Eq. (|2- lOj) . and 
using Eqs. OTT]) . (EFT!! and ([243]), we have 



P(m,n\r) = Tr ^i7 n C/(r)0- 1 0p init 77 m C/t(r; 

= Tr (n n u\T)p- mit n m u(T) 



= -^Z~ B Tr { n mU(r)n n U^T)) , (2-14) 

where we used OU(tf = U(t)G. From Eq. (12T41) . we obtain a symmetry relation 
about the joint probability as follows: 

p (m,n\T) = ^ A{ EA A _ E A A y^ E B B _ E B XB) 
P(n, m|r) 

_ Q^PQm— ingto^E m ^„ (2 - 15) 



where 



Qm^n — Em A ~ E n A ( 2 " 16 ) 



is the energy decrease in the system A, 



AE m _ n - E^ A + E% B - E mA - E^ g (2-17) 
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is the energy change of the total system caused by the on and off of the coupling 
Hamiltonian and 

A(3 = f5 B -f3 A . (2-18) 

We interpret Q m ^ n as the heat draining from the system A to the system B. Jarzyn- 
ski and Wojcik 23 ) have derived Eq. (|2-15p in the same situation. 

2.2. Deviation from the XFT and the IXFT 

Let us now retrace the derivation of the exchange fluctuation theorem 23 ) and 
show the deviation from the theorem. Using the joint probability (|2-8|) . the proba- 
bility distribution of the heat transfer Q m ^n is given by 

Pt(Q) = ^^P{ , m,n\T)5{Q - Q ), (2-19) 

where <5(-) is the delta function. Using the symmetry relation of the joint probabil- 
ity (pTo) in Eq. ([249]), we have 

PAQ) - e^V(-Q) = e Af3Q ^ P ( n , m \r) ^ AE — - l) 5 (Q - Q„) . 

m,n 

(2-20) 

Jarzynski and Wqjcik 23 ) assumed that the total energy of the systems A and 
B is almost preserved between the two measurements if the coupling 7 between the 
two systems is sufficiently weak; 

AE m ^ n = E n -E m ~0. (2-21) 

In other words, they assumed 

P(m,n\r) <x5 AEm _> nfi , (2-22) 

in the weak coupling limit. The right-hand side of Eq. (|2-20|) vanishes under this 
assumption, and thereby they obtained the exchange fluctuation theorem: 

PAQ) =i e Af3Q Pr (-Q). (2-23) 

However, we must pay attention to the behavior of the net heat transfer (Q) T 
in the weak coupling limit because we are interested in a non-equilibrium system, 
or in the presence of a finite heat transfer. From this point of view, we examine 
whether the exchange fluctuation theorem holds or not in the presence of a finite 
heat transfer for weak coupling. 

In general, the deviation term (the right-hand side of Eq. (|2-2Up ) has a finite 
value. We can check this by examining the deviation from the integral exchange 
fluctuation theorem, which is obtained by multiplying Eq. (|2-20p by and 
integrating it over the heat transfer Q as 

{e- Al3Q ) T - 1 = ^ P(n, m \ T ) (eP aAEm " n - l) , (2-24) 

771,71 
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where we used J dQ p T (Q) = 1. 

In order to compare the magnitude of the deviation from the exchange fluctu- 
ation theorem and the net heat transfer for weak coupling, we use the interaction 
picture and see the 7 dependence of the joint probability. The time-evolution oper- 
ator is written as 



U(t)=e-i Hot e- i ^ t \ 
where C(t) is a Hermitian operator defined as follows: 



(2-25) 



iC(t) 



h 



E 



t) (S Ho ) k H c . 



(2-26) 



k=o (k + iy. vr 

Here we used e A Be~ A = Y^kLo(fi-A.) k B an< ^ ^ A denotes the inner-derivation operator 

6a = [A, ] . (2-27) 

Using Eq. f)2-25|) in the transition probability (|2-9p . we have 

2 



= Tr ( n n e~^ c ^ n m e^ c ^ 



and the joint probability becomes 

P„(r) = Tr (n n e-^ c ^ Pinit n m e^ c ^ 



(2-28) 



(2-29) 



Using the joint probability in the trace form (|2-29|) . we have the deviation term from 
the integral exchange fluctuation theorem as 



(e-n - 1 = E Tr (^(W.^) (e^ AE — - 1 

m,n 

= £ Tr (e-^ZZ-^W^^I^e-^W 
-^Tr(/J m e^) Amt /7 ri ,e-^)) 



Tr p mit e 



Ho „-nC(r) p -/3 fl Ho J"fC(r) 



1. 



(2-30) 



We then expand the exponential e 1 ' 1 * 7 '^ 7 ") to obtain 

00 1 



fc=0 



i 7 Tr fpinite^ BHo fc(r),e~ feHo 
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k=2 
■ — i^Tr ^ |^ 

oo ^ 

k=2 
oo ^ 



n . ., p PbH 



(5 C{ r)) k ^ BH ° 

C(r) 



k=2 



Anite feH() (5c(r)) fc e- 



(2-31) 



where we used ^ bH ° , Pi n it] = in the last line. This commutability is a consequence 
of taking the initial state as a Gibbs state. We thus see that the deviation is of the 
second order of 7 in the lowest order. In fJH the deviation is explicitly calculated for 
simple models and is shown to have a finite value for finite 7. 

Next, we derive the net heat transfer in the power series of 7 and show that its 
lowest order of 7 is also the second order: 



Tj JlC(r)\ ( E A _ E A \ 



(Q)t = / dQ p T (Q) Q 

m,n 
m,n 

=Tr (e-^M/wtf^M - Ha^ ^ P x^° {r) 
=Tr (p-^H A - pinite^ C(r) ^e-^ C(r) 

We again expand the exponential e^* 7 ^ 1- ) to obtain 

00 1 

(Q)r =Tt Pinit H A - £ ^(^) fe Tr [ Pinit (5 c{T) ) k H A 

k=0 

00 ^ 

= - i 7 Tr (p init [C(r),H A }) - ^ ^(H) k Tr \p init (5 c{T) ) k H A 



(2-32) 



k=2 
00 



i 7 Tr (C(r)[H A , p iait \) - ^(^) fe Tr p init (5 c(T) ) k H A 



k=2 



00 _. 

^-(i7) fc Tr[p init (J c{T) ) fc ^ A ] , 



(2-33) 



k=2 



where we used [pmit,H A ] = in the last line. The lowest order of 7 in (Q) T is the 
second order, which is the same as the deviation term in Eq. (|2-3ip . This means 
that the smaller the deviation term becomes for weak coupling, the less the net heat 
transfer between the two systems becomes. We can show that the higher moments 
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of p T {Q) has the same dependence of 7, and thus the exchange fluctuation theorem 
becomes just a trivial relation in the limit 7 — > 0: p T (Q) = S(Q) = e zi ^<5(— Q) = 
5(-Q). 

In 21 we show for specific models that the deviation term (e _zi ^) T — 1 and the 
net heat transfer {Q) T are, in general, both finite in the second order of 7. 

§3. Commutable-coupling condition 

We showed in the previous section that the exchange fluctuation theorem has a 
finite deviation term in general, as far as the coupling strength 7 is finite. Before 
confirming this for specific models, we present an additional condition for which the 
exchange fluctuation theorem and the integral exchange fluctuation theorem hold 
exactly under a finite heat transfer and for arbitrary coupling strength 7. Our 
additional condition is 

[H ,H C } = Q, (3-1) 

which we refer to as the commutable-coupling condition. Jarzynski's equality for 
open quantum systems was discussed under this condition. 25 ) With the commutable- 
coupling condition, the energy of the system Hq is conserved throughout the process 
and hence AE m ^ n = holds exactly in Eq. (|2-20p . 

Let us be more precise. Under the commutable-coupling condition, we can sep- 
arate the Hamiltonian Hq and H c in the time-evolution operator, and thus the tran- 
sition probability (|2-9|) is reduced to 

T m ^ n (T) = \(n\e- i ^ I t T \m)\ 2 . (3-2) 
Using the commutable-coupling condition (|3-ip again, we have 

(E n - E m ) (n|e-^ T |m) = (n|# e~ i7 ^ T |m) - (n| e -^ r H \m) = 0, (3-3) 

and therefore we obtain 

(n\e-*r% T \m) = {u^^^ae^q. (3-4) 

We can see that the energy of the total system is preserved throughout the two-point 
measurement process. Using Eq. ()3-4p . we have 

P(m,n\T) oc T m ^ n (r) oc 5AE m ^ n ,o, (3-5) 

and thereby conclude that the exchange fluctuation theorem becomes an exact rela- 
tion under the commutable-coupling condition: 

MQ)-e A0Q pA-Q) 

= e A ^ Yl P ( n ' Ht)<W_,o (e fe ^— - l) 8 (Q - Q m ^ n ) = 0. (3-6) 

m,n 
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The integral exchange fluctuation theorem also holds exactly under the commutable- 
coupling condition, since it is a direct consequence of the exchange fluctuation the- 
orem. 

Next, we show that a finite heat transfer between the two systems does exist un- 
der the commutable-coupling condition. Under the commutable-coupling condition, 
the Hermitian operator C{t) becomes 



C(t) 



dse^ HoS H c e-^ HoS = -H c . 

n 



(3-7) 



From Eqs. (|2-33p and (|3-7p . the ensemble average of the heat transfer can be written 
in the trace form: 



(Q)r 



k=2 



k\ \ h J 



Note that 8h c Ha is generally finite, although we assume &h c Hq = 5h c {Ha+Hb) = 0. 
Therefore the heat can flow between the two systems and the exchange fluctuation 
theorem is a nontrivial relation under the commutable-coupling condition. We check 
this for specific models in the next section. 

§4. Examples 

4.1. A two-spin 1/2 system 

As the first example, we consider a quantum system which consists of two spin 
l/2s initially prepared at different temperatures. The spins exchange the energy via 
the Heisenberg coupling between the two spins. Thus, this system is given by the 
following Hamiltonian: 



H = H + 7 H C , 
H Q = H A ®l B + 1 A 



Hb, 



where Ha, Hb and H c are 



H b 
H r 



" 2 

" 2 
J ^ - 

~~^°A ■ 0~B 



(4-1) 
(4-2) 

(4-3) 



with o~ x ,a y and a z denoting the Pauli matrices. This model is analytically solvable 
and the deviation from the integral exchange fluctuation theorem, (e _zi ^) T — 1, and 
the heat transfer {Q) T are calculated explicitly as follows: 



(e 



-A/3Q\ 



2 t2 



7 V 



{e A -e B ) 2 +l 2 J 2 



sech 



'Pa^a' 


sech 


~Pb£b~ 


2 




2 
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x sinh 



x sin 



Pa) 



sinh 



Pb , 



(4-4) 



(Q)r 



£A 



2 t2 



2 (€A — £b) 2 + J 2 J 2 



sech 



Pa^A 


sech 


Pb^b 


2 




2 



x sinh 



2 (^ £ s 



x sin 2 — ^{e A -€ B ) 2 + 7 2 J 2 



2h 



(4-5) 



1 and the heat transfer (Q) T 



The r dependences of the deviation term (e 
are shown in Fig. [TJ 

Note that (e _zA/3Q ) r - 1 and (Q) T have the same dependence on 7 in this model. 
Thus, the ratio of these two quantities is independent of the coupling strength 7: 



-A/3Q 



>r"l 



sinh [<*(/J fl -/? A )] sinh 



-rye a 



(Q) T /m S mh[±(p B e B -f3Ae A )] ' (4 ' 6) 

Figure [2] shows the parameter dependence of the ratio. The ratio has a finite value 
for almost all range of the energy-level difference e B — (.a- As a consequence, if 
the coupling strength 7 is weak enough to neglect the deviation from the integral 
exchange fluctuation theorem, Eq. (|4-4|) . the net heat transfer {Q) T in Eq. (|4-5|) is 
also negligibly small. This result clearly shows that the exchange fluctuation theorem 
does not generally hold in the presence of a finite heat transfer. 

However, the ratio in Fig. [5] vanishes for €a = e_B! that is, the integral exchange 
fluctuation theorem is recovered at this particular point with a finite heat transfer. 
At this point €a = e_B and at this point only, the commutable-coupling condition is 
satisfied: 



J t 



° V A°B 



[H ,H C \ i-[<\\ - (LS 
Under the commutable-coupling condition €a 
(e -4» >T _i = o, 



°A a B 



(Q)r 



■ sech 



12 



Pa 



sech 



B 



sinh 



0. (4-7) 

e B , Eqs. flUD and ([4^5]) become 

(4-8) 

7 J. 



(Pb - Pa) 



sin 



2h 



(4-9) 



The net heat transfer has a finite value, while the integral exchange fluctuation 
theorem holds. Thus we demonstrated in this model that the integral exchange fluc- 
tuation theorem holds if and only if the commutable-coupling condition is satisfied. 

4.2. Coupled harmonic oscillators 

The second example is a system which consists of two harmonic oscillators. This 
system is given by the following Hamiltonian: 

H = H G + 1 H C , 
H = H A ®l B + 1 A 



H B . 



(4-10) 
(4-11) 
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Fig. 1. The quantities (e A/3 ®) T — 1 (solid red line) and (Q)t (broken blue line) in the two-spin 
1/2 system. We fixed p A e A = 2: for e B /e A = 0.8, Pb/Pa = 0.9 with (a) -yJ/e A = 0.1 and (b) 
jJ/e A = 0.01; for e B /e A = 0.8, Pb/Pa = 1.1 with (c) -yJ/e A = 0.1 and (d) -yJ/e A = 0.01; for 
e B /e A = 1.25, p B / Pa = 0.9 with (e) 1 J/e A = 0.1 and (f) 1 J/e A = 0.01; for e B /e A = 1.25, 
Pb/Pa = 1.1 with (g) <yJ/e A = 0.1 and (h) <yJ/e A = 0.01. 



Here the Hamiltonian Ha, Hb and H c are 

Hb = foujB^b, 
H c = v [a^b + tfa 

where a and a) are the creation and annihilation operators of the oscillator A, b and 
&t are those of B, [a, a'] 
number. 



(4-12) 



1, [a, a] = and [b,bM = 1, [b,b] = 0, and v is a real 
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Fig. 2. The ratio of the (e~ Al3Q ) T - 1 to (Q) r . We fixed /3 A e A = 2 with (a) /3 s //3a = 0.9 and (b) 
(3b /Pa — 1.1. The commutable-coupling condition is fulfilled at the point es — £A = 0. 



To the second order of the coupling strength 7, the deviation from the inte- 
gral exchange fluctuation theorem, [e~ A ^) T — 1, and the heat transfer (Q) T are 
analytically calculated as follows: 



(e- 4 «) 



1 =4 



ju sin [|(aM - uj b )} 1 
ft (wi - cjb)/2 J 

sinh (Af3hu>A/2) sinh [/^^(cja — <^_b)/2] 



2 sinh {Mp. 



2sinh(&Mp) 



+ 0( 7 3 



(4-13) 



71/ sin 



I ft (wa - wb)/2 J 

^ sinh - /? j b^u; s )/2] 



2 sinh 



2 sinh 



+ 0( 7 3 



(4-14) 



The r dependence of the deviation from the integral exchange fluctuation theorem 
and the net heat transfer during r are shown in Fig. O To the second order of the 
coupling strength 7, the ratio of the above two quantities is given by 



sinh (A[3Huja/2) sinh \J3b^{<^a — ^b)/2] 
sinh [h(f3 B UB ~ (3aua)/2\ 



+ 0( 7 3 



(4-15) 
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Fig. 3. The quantities (e ) T — 1 (solid red line) and (Q} T (broken blue line) in the coupled 
harmonic oscillators. We fixed PaHuja = 4: for ojb/ua = 0.8, (3b /Pa = 0.9 with (a) jv/Huja = 
0.1 and (b) -)v/hu A = 0.01; for lub/uja = 0.8, /3b//3a = 1.1 with (c) -yu/hu) A = 0.1 and (d) 
-yv/tkuA = 0.01; for lub/loa = 1.25, fa /Pa = 0.9 with (e) -jv/huA = 0.1 and (f) -yv/htUA = 0.01; 
for ljb/ua = 1.25, Pb/Pa = 1.1 with (g) ^u/Kuja = 0.1 and (h) 'yv/fkJA = 0.01. 



This result again shows that the exchange fluctuation theorem does not generally 
hold in the presence of a finite heat transfer. 

However, we can see that the integral fluctuation theorem holds in the presence 
of a finite heat transfer at the point given by uja = u>b- At this point and at this 
point only, the commutable-coupling condition is satisfied as 

[H , H c ] = hu{to A - uj B )(a^b - 6 f a) = 0. (4-16) 

For the more general coupling Hamiltonian, H c = v [(a^) lA b lB + (b') lB a lA ], the 
deviation from the integral exchange fluctuation theorem and the net heat transfer 
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are calculated to the second order of 7 as 

(e -A(3 Q) 1 _ $ 71s sin [^(uaIa-ubIb)] V 
{ h \h {oo a Ia-ujbIb)/2 J 

4/ A ! Ib\ 



[2smh((3 A hio A /2)} lA [2sinh(/? B ^; B /2)]' s 
x sinh {Af3hio A lA/2) sinh [/3 B h(L0 A l A - u B l B )/2] + 0( 7 3 ), (4-17) 



(Q)r = 



sin [f (u^A ~ u b Ib)] } 
h (lj a Ia- uj b I b )/2 J 

-2l A \ Ib\ 



2 



[2 sinh (/?a^a/2)]^ [2 sinh {(3 B hw B /2)} 1 



x IUoaIa sinh [(/a/?a^a - lBpBhu B )/2] + 0( 7 3 ). (4-18) 
The ratio of the above two quantities is given by 

(e~ A P Q ) T - 1 _ sinh (Af3huj A lA/2) sinh [(3 B h{oj A lA -uj b I b )/2] 



{Q)r/{^4 A ) sinh WbwbIb ~ PaIaua)/2] 



+ 0( 7 3 ). (4-19) 



This shows that the integral exchange fluctuation theorem holds in the presence of 
a finite heat transfer if and only if the commutable-coupling condition is satisfied: 

[H , H c ] = hu{uo A l A ~ w B I B ) [{o)) lA b lB - {tf ) lB a lA \ = 0, (4-20) 

or luaIa = u B l B . 

§5. Conclusions 

To summarize, we showed that the exchange fluctuation theorem in its original 
form does not generally hold in the presence of a finite heat transfer. In the limit 
7 — ► 0, the kth. moments of p T (Q) vanish for k > 1. The deviation from the exchange 
fluctuation theorem also vanishes in the same order of 7. We derived general formulas 
for the above and analytically confirmed them for specific models. This means that 
there is no heat transfer when the coupling strength 7 is small enough to neglect 
the deviation from the exchange fluctuation theorem. In this case, the exchange 
fluctuation theorem reduces to a trivial relation and has no information about the 
heat transfer. 

However, we found a condition for the exchange fluctuation theorem to hold 
exactly for arbitrary 7 and we referred to it as the commutable-coupling condition. 
Under this condition, the exchange fluctuation theorem becomes an exact relation 
independently of the coupling strength 7 under the existence of a finite heat transfer. 
We confirmed this in specific models. In short, the exchange fluctuation theorem 
holds not because the coupling is weak as was originally proposed, but because the 
total energy of the system is conserved under the commutable-coupling condition. 



16 



S. Akagawa and N. Hatano 



The deviation from the exchange fluctuation theorem consists of the commuta- 
tion relation between the Hamiltonian of the total system and the coupling Hamil- 
tonian. Therefore, the non-commutativity of the observable in quantum mechanics 
plays an important role in the deviation. 
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